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QCDSF/UKQCD COLLABORATIONS 

Generalized transversity distributions encode essential information on the internal 
structure of hadrons related to transversely polarized quarks. Lattice QCD allows 
us to compute the lowest moments of these tensor generalized parton distributions. 
In this talk, we discuss a first lattice study of the Soffer bound and show preliminary 
results for transverse spin densities of quarks in the nucleon. 



1. Introduction 

Generalized parton distributions (GPDs) [1,2] are an ideal tool to study 
many fundamental facets of hadron structure in terms of quarks and glu- 
ons. One key point is the relation of GPDs to (orbital) angular momentum, 
which plays a central role for the nucleon spin sum rule [3] . Moreover, GPDs 
allow us to investigate the nontrivial interplay of longitudinal momentum 
and transverse coordinate space degrees of freedom [4-6]. In this contri- 
bution, we will focus our attention on the recently observed correlation of 
transverse quark spin and impact parameter which shows up in transverse 
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spin densities of quarks in the nucleon [7]. It turns out that these correla- 
tions in the transverse plane are governed by quark helicity flip (or tensor) 
GPDs [8]. As in the case of the unpolarized and the polarized GPDs [2], 
they are defined via the parametrization of an off-forward nucleon matrix 
element of a bilocal quark operator as follows 

(P',A'\ j°° ^e iXx q(-±r?) n^ l5 Uq(±n^ \P, A) = u(P\ A>„{ 
^ 75 [h t {x, £ t) ^_H T (x, *)) + t^OS. E T {x, £, t) 

+ 2 J Ht{x > * } + m 7 " ^(x, t)J tt (J> A) , (1) 

where /I'H = ft 1 " — A = P' — P is the momentum transfer with 
t = A 2 , P = (P' +P)/2, and £ = — n ■ A/2 defines the longitudinal momen- 
tum transfer with the light-like vector n. The Wilson line ensuring gauge 
invariance of the bilocal operator is denoted by It. Our parametrization 
in Eq. (1) in terms of the four independent tensor GPDs slightly differs 
from the literature [7,8] where a function Et instead of Et has been used. 
However, in [7, 9] it has been noted that Et typically appears in linear 
combination with the tensor GPD Ht- It is therefore reasonable to adopt 
a new notation and consider Et = Et + 2Ht and Ht as fundamental 
quantities. 

One prominent feature of GPDs is that they reproduce the well known 
parton distributions in the forward limit, A = 0, and that their integral 
over the momentum fraction x leads directly to form factors. For this rea- 
son, the GPD Ht(x, £, t) is called generalized transversity, since for vanish- 
ing momentum transfer it is equal to the transversity parton distribution, 
H T (x,0,0) = Sq(x) = fn(x) for x > and H T (x,0,0) = -Sq(-x) = 
—hi(—x) for x < 0. On the other hand, integrating Ht{x, £, t) over x gives 
the tensor form factor: 

J dxH T {x,£,t)=g T {t). (2) 

Another feature of GPDs important for our investigations below is their 
interpretation as densities in the transverse plane for £ = [4]. To give an 
example, it has been shown that the impact parameter dependent quark 
distribution for the quark GPD H qi 

q(x,b ± ) = J ^e- ib -- A -H q (x,Z = 0,t = -Al), (3) 
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has the interpretation of a probability density for unpolarized quarks of 
flavor q with longitudinal momentum fraction x and transverse position 
b± = (b x , by) relative to the center of momentum in a nuclcon. 

In order to facilitate the computation of the tensor GPDs in lattice 
QCD, we first transform the LHS of Eq. (1) to Mellin space by forming the 
integral j_ x dxx n ~ 1 ■ ■ ■ . This results in nucleon matrix elements of towers 
of local tensor operators 

0^i-^-i(O) = q(0)ia^ l5 iD^ . . . S^-^q(O) , (4) 
which are parametrized in terms of tensor generalized form factors (GFFs) 

Ami, B Tn i, A Tn i and B Tn i- Here, D = ^(D — D) and {•••} indicates 
symmetrization of indices and subtraction of traces. a For n = 1, we have 
[8,10] 

(P'A'| 9(0)^755(0) |PA) = u(P',A')^ v 75 (A TW (t) - J^A TW {t)j 

+ 2 t + 2 J ] A ™(t)\u(P,A) . (5) 

The relation of the lowest moment of the tensor GPDs to the GFFs is 
simple and given by 

HfTH&t) =A TW {t) = g T (t), HT\Z,t) = A TW {t) 
E^&t) = B TW (t), E™ =1 (^t) = B TW {t) =0, 

where iJ^(^,t) = J^ 1 dxx n ~ 1 HT(x,l;,t). The general parametrization in 
terms of GFFs and their relations to the moments of the GPDs for n > 1 
can be found in [10, 11]. 

The calculation of moments of GPDs in lattice QCD follows standard 
methods, which have been described in detail in the literature [12-14]. In 
the following, we therefore give only an outline of the procedure we use 
to extract the GFFs. First, nucleon matrix elements in the form of two- 
and three-point functions are computed on an Euclidean space-time lat- 
tice. The typical suppression of the matrix elements by exponential factors 
exp(— tE) in the Euclidean time r and the energy E is cancelled out by con- 
structing an appropriate ratio R{t) of three- to two-point functions, which 
is averaged over the plateau-region i?(r p i a t.) ~ const. The averaged ratio 
is then renormalized and equated with the continuum parametrization of 
the corresponding nucleon matrix element, e.g. Eq. (5), for all contributing 



a The Mellin-moment index n used in this work differs from the n in [10] by one. 
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Figure 1. Lattice results for the lowest two moments of the Soffer bound for up-quarks. 

index (/i, v) and momentum (P, P') combinations. This leads to an overde- 
termined set of linear equations which is solved to extract the GFFs. The 
statistical error on the GFFs is obtained from a jackknife analysis. Our 
results have been non-perturbative renormalized [15] and transformed to 
the MS scheme at a scale of 4 GeV 2 . 

The lattice results to be discussed below have been obtained from simu- 
lations with rif = 2 flavors of dynamical non-perturbatively 0(a) improved 
Wilson fermions and Wilson glue. There are 12 datasets available consist- 
ing of four different couplings (3 — 5.20, 5.25, 5.29, 5.40 with three different 
k = «soa values per j3. The pion masses of our calculation vary from 550 to 
1000 MeV, and the lattice spacings and spatial volumes vary between 0.07- 
0.11 fm and (1.4-2.0 fm) 3 respectively. Our calculation does not include 
the computationally demanding disconnected contributions. We expect, 
however, that they are small for the tensor GFFs [14]. More details of the 
simulation can be found in [14, 16, 17]. 



2. Lattice study of the Soffer bound 

In this section, we investigate the Soffer bound [18] 

\6q(x)\ <±(q(x) + Aq(x)) , (6) 

which holds exactly only for quark and anti-quark distributions separately. 
For discussion of its validity, see e.g. [19] and section 3.12.3 of [2]. For a 
lattice study of the Soffer bound, we take Mellin moments of Eq. (6) and 
consider the "Soffer-ratio" 



S n = 



2|(x"- 1 )^| 
(x™" 1 } + {x n - 1 ) A 



n=l,2... 



(7) 
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where (x rl_1 ) = J_ x dxx n ~ 1 q(x). At this point it is important to note that 
Mcllin moments of distribution functions give always sums/differences of 
moments of quark and anti-quark distributions, e.g. = (x n ^ 1 ) q + 

(—l) n {x n ^ 1 )q. Therefore, the ratio S n in Eq. (7) is not necessarily smaller 
than one. Experience shows that contributions from anti-quarks are neg- 
ligible in our calculation. In Fig. (1) we show our results for the ratio (7) 
versus the pion mass for up-quarks (similar results for down-quarks can be 
found in [14]). The fact that the ratio is consistently below one for the low- 
est two moments of the up and the down quarks strongly suggests that the 
Softer bound is satisfied in our lattice calculation. The lattice results show 
almost no dependence on the pion mass due to cancellations of the pion 
mass dependence of the individual distribution functions in the ratio (7). 
Linear chiral extrapolation in m 2 . leads to the following predictions for the 
ratios at the physical pion mass 

up-quarks : 5" =1 = 0.60±.01, S n = 2 = 0.78±.01 
down-quarks: 5" =1 = 0.57±.02, S n = 2 = 0.73±.05 . " M 



3. Lattice results for the lowest moment of the transverse 
spin density 

We now turn our attention to a discussion of our lattice results for the 
density of transversely polarized quarks in the nucleon. The lowest moment 
of the quark transverse spin density is given by [7] 

(P + ,R± = 0,S±\ \q{bi_)[ 1 + - s j ± ia +j ^ 5 ]q(b ± ) |P+ R ± = 0, S±) = 

\ |ai (6_l) + s\S\ (A T1Q (b ± ) - JL^A b± A TW (b±)j 

+ b ^(siB[ (b ± )+s i 1 B' TW (b ± j) 

+ s i(2&i&i - bi^si^A''^)} . (9) 

The transversity states 

|P+ R ± = 0, S±) = ( \P + , R± = 0, A = +)+e 4X |P+ R ± = 0, A = -) ) 

2 (10) 
describe a nucleon with longitudinal momentum P + = (P° + P 3 )/v / 2 which 
is localized in the transverse plane at R± = and has transverse spin 
Sj_ = (cos x, sin x). The impact parameter dependent GFFs in Eq. (9) are 
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Figure 2. Densities of up-quarks in the nuclcon. The nucleon and quark spins are 
oriented in the transverse plane as indicated, where the inner arrow represent the quark 
and the outer arrow the nucleon spin. A missing arrow represents the unpolarized case. 



just the Fourier-transforms of the momentum space GFFs at £ = 0, as in 
Eq. (3). The derivatives in Eq. (9) are defined by f'(b±) = d b 2 f(b±) and 
A b± f(b ± ) = id bl (bld bl )f(b ± ). 

Since momenta are discretized on a finite lattice, we obtain the GFFs 
only for a limited number of different values of the momentum transfer 
squared t. We have in general 16 i-values available per dataset in a range 
of < t < 4 GeV 2 . To facilitate the Fourier transformation to impact 
parameter space, we parametrize the GFFs using a p-pole ansatz 

(1 - t/ml) 

where the parameters F(0), m p and p for the individual GFFs are fixed by 
a fit to the lattice results. The ansatz in Eq. (11) is then Fourier trans- 
formed in order to get the GFFs as functions of the impact parameter b±. 
Details of the p-pole parametrization and numerical results for the param- 
eters will be given in a separate publication [20]. Here we only note that 
the values for the power p we are using for the different GFFs lead to a 
regular behavior of the transverse spin density in the limit b± — > 0, as 
discussed in [7]. In Figs. (2,3) we show preliminary results for the lowest 
moment of transverse spin densities of quarks in the nucleon for up quarks. 
We note that the plots do not exactly show probability densities because 
the lowest moment corresponds to the difference of quark and anti-quark 
densities. The densities are however strictly positive for all b±, indicating 
that the contributions from anti-quarks are small. On the LHS of Fig. (2), 
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Figure 3. Up-quark densities. Symbols are explained in the caption of Fig.(2). 

we show the transversely distorted density of unpolarized quarks in a nu- 
cleoli with spin in x-direction (coming from the dipolc-tcrm oc e Jl hP^_S\ in 
Eq. (9)), which has already been discussed in [21]. A new observation is 
that the GPD Et also leads to a strong transverse distortion orthogonal 
to the transverse quark spin for an unpolarized nucleon (coming from the 
dipolc-tcrm oc e> l \?^s\ in Eq. (9)) on the RHS of Fig. (2). It was argued 
in [9] that this shift in +j/-direction b may correspond to a non-zero, nega- 
tive Boer-Mulders function [22] < for up-quarks. The distortions due 
to transverse quark and nucleon spin add up for the density on the LHS 
in Fig. (3), while it goes in opposite (— j/)-direction for quarks with spin 
opposite to the nucleon spin, as can be seen on the RHS of Fig. (3). Inter- 
estingly, there is practically no influence visible from the quadrupolc-term 
oc s\(26\6^ - b'j_5 lj )S J L in Eq. (9) for the up-quark densities. 

4. Conclusions and outlook 

Our lattice results for the transversity distribution suggests that the Soffer 
bound is saturated by w 60 — 80% for the lowest two rr-moments. In ad- 
dition, we have presented preliminary results for the lowest moment of the 
transverse spin density of quarks in the nucleon. The distortion of the den- 
sity of transversely polarized quarks in an unpolarized nucleon is substantial 
and could give rise to a non-vanishing negative Boer-Mulders function for 
up-quarks through final state interactions as argued by Burkardt [9]. 



b or equivalently a shift in (— ^-direction for quarks with spin in y-direction 
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We plan to extend our analysis of transverse spin densities in lattice 
QCD to the lowest two moments of up- and down-quarks and to investigate 
improved positivity bounds for GPDs which have been obtained in [7] . 
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